TO APPEAR IN IEEE TRANSACTIONS ON SIGNAL PROCESSING 1 

An Analytical Framework for Heterogeneous 
Partial Feedback Design in Heterogeneous 
Multicell OFDMA Networks 

Yichao Huang, Member, IEEE, and Bhaskar D. Rao, Fellow, IEEE 

Abstract 

The inherent heterogeneous structure resulting from user densities and large scale channel effects 
motivates heterogeneous partial feedback design in heterogeneous networks. In such emerging networks, 
a distributed scheduling policy which enjoys multiuser diversity as well as maintains fairness among users 
is favored for individual user rate enhancement and guarantees. For a system employing the cumulative 
distribution function based scheduling, which satisfies the two above mentioned desired features, we 
develop an analytical framework to investigate heterogeneous partial feedback in a general OFDMA- 
based heterogeneous multicell employing the best-M partial feedback strategy. Exact sum rate analysis 
is first carried out and closed form expressions are obtained by a novel decomposition of the probability 
density function of the selected user's signal-to-interference-plus-noise ratio. To draw further insight, we 
perform asymptotic analysis using extreme value theory to examine the effect of partial feedback on 
the randomness of multiuser diversity, show the asymptotic optimality of best- 1 feedback, and derive an 
asymptotic approximation for the sum rate in order to determine the minimum required partial feedback. 
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I. Introduction 

The growing dependence of users on wireless services will require wireless systems to become ubiqui- 
tous and offer seamless support. The demands of video and other high data rate applications have placed 
increasing requirements on networks to support high data rate services in a cost effective manner leading 
to heterogeneous networks. With the advent of OFDMA-based heterogeneous networks fll], El which 
incorporate lower power pico ||3], [@], femto base stations [0, 0, and fixed relays Q, (H to coexist 
with the traditional macrocell, the spectrum reuse has grown aggressively to full usage pattern across 
different tiers [9 1 of the heterogeneous structure. One challenging feature of heterogeneous networks is 
the self-created "cell edges" within the macrocell, which requires advanced techniques both in theory 
and in practice to model |10], manage ifTTI . and even make use of the cross-tier intercell interference. 
Among the recent approaches that have been developed such as subcarrier allocation and power control to 
carefully adapt the system resource in a centralized way lli2l - |[T4ll , and utilization of spatial domain for 
cooperative multicell processing to cancel, coordinate, and align the interference |fT31 -|20|, a significant 
impediment is the NP-hardness of the problem ET1 . the limited resource constraints, as well as the 
need for extensive backhaul capability. These challenges favor the development of distributed solutions. 
All the aforementioned techniques in the downlink assume the availability of channel state information 
at the transmitter (CSIT) via feedback [22] to adapt the network transmission strategy to the varying 
wireless environment. With the rapidly growing number of wireless users, the amount of feedback for 
the OFDMA-based networks may become prohibitive which motivates the design of efficient feedback 
schemes without significantly degrading system performance. 

In addition to the usual challenges, there are two new issues that arise when investigating partial 
feedback in heterogeneous networks. Firstly, due to the different locations of the users and different 
ranges of transmit powers, users' large scale channel effects are highly asymmetric. Therefore, it is 
equally important to guarantee fairness among users as well as leveraging multiuser diversity [23], Il24l 
in an opportunistic scheduling framework. Secondly, since the number of users or user densities are diverse 
in a heterogeneous network, it would be beneficial to adapt the feedback and require less feedback when 



'For the purpose of scheduling or optimization in a multicell network, the feedback is often needed in a frequency division 
duplex (FDD) system, or in a time division duplex (TDD) system when the channel reciprocity can not be observed due to 
duplexing time delay. 
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a serving base station has more users associated with it. We refer to this methodology as heterogeneous 
partial feedback and we aim to provide an analytical framework to identify its benefits under a fair and 
distributed opportunistic scheduling policy to fulfill the vision of location awareness E71 and situational 
awareness. 

The first step towards examining the aforementioned issues is developing an opportunistic scheduling 
policy, which exploits multiuser diversity and also preserves scheduling fairness among heterogeneous 
users. Traditional scheduling policies such as round robin strategy [28 ] and greedy strategy [23 ] are easy 
to implement, yet only achieve one of the desired features. In this paper, we consider a system that 
employs the cumulative distribution function (CDF) based scheduling policy 1291 , [30]. According to the 
basic CDF-based scheduling strategy, an user is selected whose rate is high enough, but least probable 
to grow higher. Therefore, this scheduling strategy possesses properties similar to the proportional fair 
scheduler ll24l . ||3TI - ||341 . and additionally enables a user's rate to be independent of the statistics of 
other users. Herein, the CDF-based scheduling policy is analyzed for a general OFDMA downlink in a 
multicell environment to examine heterogeneous partial feedback design. Currently, in an OFDMA-based 
system which groups subcarriers into resource blocks l55l - lf37l to form the basic scheduling and feedback 
unit, two partial feedback strategies are appealing: the thresholding-based partial feedback [38 1 — [40 1 and 
the best-M partial feedback PTI - ll46l . The latter strategy, which is considered in practical systems such 
as LTE |[25l . [26], requires the users to order and convey the M best channels. Herein, we employ the 
best-M partial feedback strategy for further analysis. Intuitively, M would be chosen to be small when the 
user density in a given cell is large, which motivates the utilization of heterogeneous feedback resource 
across different cells in the heterogeneous networks. 

Rigorous development of the analytical framework requires investigation of the interplay between the 
scheduling policy, partial feedback, and the statistical property of the user's signal-to-interference-plus- 
noise ratio (SINR). There are limited results available in the literature on this topic and the available 
results analyzing the best-M partial feedback are for the single cell scenario without intercell interference 
BT1 -I46I. A detailed treatment of the best-M partial feedback is provided in [46] and a convenient 
polynomial form for the CDF of selected user's SINR is presented for analytical evaluation. Though it is 
derived for the single cell scenario, it forms the building block for our heterogeneous network analysis 
which takes into account the cross-tier intercell interference. In this paper, the analytical framework is 



2 Note that achieving adaptive feedback according to the number of users, the users' channel condition, and users' data rate 
requirements etc is an important issue in practical systems such as LTE [25 1, |26|. 
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treated first from the perspective of exact system performance. We derive the closed form expression 
for the sum rate with the CDF-based scheduling policy and best-M partial feedback strategy. One key 
technique developed and utilized is the decomposition of the probability density function (PDF) of the 
selected user's SINR, which is amenable for further integration needed to determine system performance. 
The derived closed form results are directly applicable to further system evaluation. 

In order to gain additional insight, we investigate the system performance from the asymptotic perspec- 
tive utilizing extreme value theory [47 1, [48 1 when the number of users in a given cell grows large [49 1- 
|54|| . Different from the special case of full feedback, examining the general best-M partial feedback incurs 
additional difficulties due to the two-stage maximization resulting from partial feedback and scheduling 
policy, with the first stage maximization being performed at the user side to provide selective feedback 
and the second stage maximization being performed at the scheduler side for user selection. Herein, we 
analyze the tail behavior of the selected user's SINR and establish the type of convergence in order to 
examine the effect of partial feedback on the randomness of multiuser diversity, show the optimality of 
best-1 feedback in the asymptotic sense, and more importantly, provide the asymptotic approximation 
for the sum rate of the general best-M partial feedback. The established asymptotic results further help 
in analytically tracking and determining the minimum required partial feedback. 

To summarize, the contributions of this paper are threefold: a conceptual framework for situational - 
aware heterogeneous partial feedback design in an OFDMA-based heterogeneous multicell network, a 
thorough analysis and derivation of closed form results for the sum rate, and a detailed investigation of 
the partial feedback based on extreme value theory. All these contributions foster the understanding of 
heterogeneous feedback design in future systems. Furthermore, the analytical tools developed promise 
to have broad applicability and can be applied to many related problems. The remainder of the paper is 
organized as follows. The system model is provided in Section ITT] The general treatment without specific 
channel models is examined in Section UlTl By assuming standard channel models, Section HVl carries out 
exact performance analysis, and Section |V] presents asymptotic analysis. Numerical results are provided 
in Section |VlJ and Section IVIII concludes the paper. 

II. System Model 

We consider the downlink of an OFDMA-based heterogeneous network. The model assumed is generic 
and sufficiently general to be applicable to a multitier multicell network e.g., see Fig. Q] for illustration. 

3 The special case with one picocell inside a macrocell under symmetric large scale channel effects is studied in our recent 
work (55). 
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Fig. 1. Illustration of a generic OFDMA-based multicell heterogeneous cellular networks. Each cell has users associated with 
it. One selected user for transmission in one resource block is shown: solid line (desired signal); dashed line (potential intercell 
interference). 



The system consists of N resource blocks, with one resource block as the basic feedback and scheduling 
unit. Full spectrum reuse is assumed and it is also assumed that there is no advanced technique employed 
to suppress interference such as multiuser detection at the receiver side. The process of cell association 
is assumed to be performed in advance. Without loss of generality, one base station Bq from the base 
station set B and its associated users JCq = {1, . . . , k, . . . , Kq} with |/Co| = Kq are considered. 
The received signal of user k at resource block n is represented by 



(o) 



(1) 



where the superscript indicates the base station, i.e., the serving cell and the interfering cells; denotes 
the number of effective interfering cells for user k, with the influence of other interfering cells, namely 
the residual interference, included in the additive white noise v^ n distributed with CJ\f(0,a^,). s^ n and 



are the transmitted symbols by the serving cell and the interfering cell B^ with E 



and E 



M |2 



,(0) |2 



P 



(0) 



P 



Hfcn anc ^ ^k ni which are assumed to be independent across users and resource 
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block^j, denote the small scale frequency domain channel transfer function between the serving cell 
and user k at resource block n, and between the interfering cell and user k at resource block n, 
respectively. GjP and G$ represent the large scale channel gain between the serving cell and user k, 
and between the interfering cell B& and user k respectively. Based on the aforementioned assumption, 

the SINR of user k at resource block n can be written as 

r(°)«(0)irr(°)|2 _(0W(0)|2 

C|Mr(0) _ U fc P \ n k,n\ _ Pk \ n k,n\ (r> . 



2^=1 P ( >\ H k,n\ +(T k 2^b=lPk \ H k,n\ + 1 



where = Gk % - - , = Gfc ^ - . The SINR is the channel quality information (CQI) that will be fed 
back and used for scheduling as discussed next. 

III. General Sum Rate Analysis with CDF Scheduling Policy and 

Best-M Partial Feedback 

This section is devoted to the analysis of the interplay between the scheduling policy and partial 
feedback for a general channel model (note that no assumption on the distribution of the large and small 
scale channel gains has been made so far), with treatment of specific channel model in Section [TV] 

Let Z^' represent SINR^ for notational simplicity and denote it as the CQI of user k at resource 
block n with CDF F z m. The CDF does not depend on the resource block index n because the SINR^'s 
are independent and identically distributed (i.i.d.) across resource blocks n for a given user k. Now 
consider the feedback procedure utilizing the best-M partial feedback strategy. According to the best-M 
partial feedback strategy, users measure CQI for each resource block at their receiver and feed back the 
CQI values of the M best resource blocks chosen from the total of N valuejfl More details on the best-M 
partial feedback approach can be found in B31 - E61 . ll56l . This selective feedback procedure involves 
a maximization stage at each user. Because only a subset of the ordered CQI are fed back, from the 
perspective of the scheduler (i.e., the serving base station Bq), if it receives feedback on a certain resource 
block from a user, it is likely to be any one of the CQI from the ordered subset. We now aim to find the 
CDF of the CQI seen at the scheduler side as a consequence of partial feedback. Denote m as the 
received CQI at the scheduler for user k at resource block n under best-M partial feedback, which is the 
outcome of the user side maximization. Let F v w be its CDF, with resource block index n dropped due 



4 This assumption corresponds to the frequency domain block fading channel model |36|, |45|, |46| due to its simplicity and 
capability to provide a good approximation to actual physical channels. 
5 We assume the CQI is fed back without feedback delay. 
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to the i.i.d. property across resource blocks for a given user. It is easy to see that for the full feedback 
case, i.e., M = N, F Y w = F z (<>), and for the best-1 feedback case, i.e., M = 1, F Y w = (F Z ( 0) ) N . 



1 k,N 



Utilizing the results in ||46l , the CDF for the general best-M feedback case can be expressed as 

M-l 

F Y m (x) = V UN,M,m)(F zm (x)) N - m , (3) 

I h,M ^— ' 

where &(N,M,m) = Y%L™ O ("I)'"™ 

After the scheduler receives feedback from its serving users, it is ready to perform scheduling. It is 
clear that for the single cell scenario without intercell interference, the scheduling policy is easier to 
implement and analyze. For instance, in the single cell scenario with homogeneous users, namely same 
large scale effects, the greedy scheduler or the max-SNR scheduler makes full use of multiuser diversity 
as well as guarantees fairness due to the same statistics of the user's CQI. In the single cell scenario 
with heterogeneous users, i.e., different large scale effects l46l . the normalized greedy scheduler which 
selects user according to their normalized CQI has the same desired property. However, in the general 
multicell scenario with intercell interference, the SINR^'s are independent but non-identically distributed 
(i.n.i.d.) across users. Therefore, the received CQI at the scheduler for different users Y^^ m i.n.i.d. 



across users. In order to leverage multiuser diversity and guarantee fairness^, we employ the CDF-based 
scheduling policy |29l , ll30l . According to the CDF-based scheduling policy, the scheduler will utilize 



the distribution of the received CQI for each user, i.e., F v (o) . Herein, it is assumed that the scheduler 
perfectly knows the CDIrj, and it conducts the following transformation 

Y k%M = F Y<% ( Y k°l,M)- ( 4 ) 

The transformed random variable M is uniformly distributed over the range from to 1 and can be 
regarded as the virtual received CQI of user k at resource block n. The transformed random variables 
Y k°l m' s are i-i-d- across users, which enables the maximization at the scheduler side to perform fair 

6 Note that the motivations as well as the fairness for the proportional-fair (PF) scheduling policy and CDF-based scheduling 
policy are very different. The PF policy targets the system utility as the definition of system fairness. The CDF-based policy 
targets the long-term user fairness and each user on average is equiprobable to be scheduled. In the single cell case, it can be 
shown that these two scheduling policies have similar effects. However, in the general multicell case, users' rates are coupled 
under the PF policy, and independent under the CDF-based policy. Analyzing and comparing these two scheduling policies in 
the general multicell networks are left to our future work. 

7 This is the only system requirement to perform CDF-based scheduling, and the CDF can be obtained by infrequent feedback 
from users and learned by the system. Methods to estimate the CDF can be found in |29|. 
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scheduling. Denoting k* n as the random variable representing the selected user for transmission at resource 
block n, then 



k* =arg max 

kEU„,M ' ■ 



(5) 



where U Hi m denotes the set of users who convey feedback for resource block n. It can be easily seen 
that when M = N, P(|W n ,Ar| = Ko) = 1. For the general case when 1 < M < N, the probability mass 
function (PMF) of \U n ^M\ can be shown to be 



M\ T ° f M 
iV j [ ~ 77 



K —T 



,0<t <K . 



(6) 



After the user k* n is selected according to d5j, the scheduler utilizes the corresponding Y^) n M for rate 
matching of the selected user. We denote the random variable X^ M as the selected user's CQI for 
resource block n, and use the sum rate as the system performance metric. The sum rate C^\M) for a 
given base station Bo employing the CDF-based scheduling and best-M partial feedback is defined as 
follows 



1 N 

n=l 



(0) 
M 



(7) 



From the aforementioned analysis, the sum rate can be formulated, with appropriate conditioning, a 



1 



N 



'|W„, M | 



E 



x . 



log 2 ( 1 + X% 



\U n , M \ + 



+ Ey(H] 



log 2 1 + X 



(0) 
n,M 



\Un, 



n=l 
N 



(a) 



/ log 2 (l + t)d(F yC o, (t)) To | \U M \^0 
Jo k *- M 



(8) 



where (a) follows from the identical distributed property across resource blocks and the change of variable 
x = F v (o) (t). The conditional statistical property of X^\ /f conditioned on the selected user k* and the 
set of users who have conveyed feedback U n ^M can be expressed as 



X^\k*=k,\U M \=T 



(x) = (F Y (o, (x)) 



Using ([3]), it can be expressed in the following power series expansion |46l , 11571 

To(M-l) 

F x£)|fc , =fe>|WM|=To (x)= UN,M,r ,m)(F zi0) (x)) N ^- m , 

m=0 



(9) 



(10) 



8 In order to maintain full frequency reuse for analytical tractability, it is assumed that if no user provides CQI for a certain 
resource block, then that resource block would be in outage and would not contribute to the sum rate calculation. 
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where 



£ 2 (N, M,r ,m) 



(£i(iY,M,0)) ro , m = 



E minim, M— I),, -. s „ \ 

l=i ; ((r + l)£-m) 



mTi (N,M,0) 

xCi(N,M,£)^ 2 (N,M,T ,m- £), 1 < m < r (M - 1) 
fa(N,M,M - 1)) T0 , m = T (M-l). 
Using © and (TTOb . the sum rate © can be expressed in the following form 



(11) 



C*°)(M) =E fc .E, WM | 



(a) J_ 
^0 



EE 

fc = l T =l 



log 2 (l +a?)dFy(o )|fc , =fc) | Wm | =T0 (x) I \U M \ + 

i 

E C2(iV,M,T 0) m)g fe (iVTo 



#0 
TO 



M 

iV 



771 



(12) 



m=0 



where (a) follows from the fair property of the CDF-based scheduling policy: P (k* = k, \Um\ = to) 
-^P (\Um\ = tq). The integration £7&(e) for e G N + is defined as 



Qk{e) 



\og 2 (l + x)d(F zm (x)Y. 



(13) 



From (fT2l . the individual user rate for user A; can be expressed as 

K -r ro(Af-l) 



= ^ E 

T = l 



i^o 

TO 



MV" / M 



N J 



N 



E 6(iV,M,To,m)g fe (iVTo-m). (14) 



m=0 



For the special full feedback case, the sum rate becomes C(°\N) = Y^k=l @k(Ko)> an d the individual 
user rate for user k becomes C^\n) = -g-Gk(Ko). 

Remark: A few remarks are in order. Firstly, the effect of best-M partial feedback and the CDF- 
scheduling policy result in a two stage maximization. The first stage maximization occurs at each user 
side to select the M best CQI for feedback. The second stage maximization is conducted at the scheduler 
side by performing CDF-based transformation and user scheduling. Secondly, with the help of CDF-based 
scheduling, each user feels as if the other users had the same CDF for scheduling competition |[29l . In 
other words, each individual user's rate is independent of other users. This important feature not only 
enables the distributed system to enjoy multiuser diversity, but also makes it possible to consider or 
predict each user's rate by only considering its own CDF. Thirdly, users are equiprobable to be scheduled 
despite of their heterogeneous channels (e.g., different statistics due to diverse propagation environments 
and interference levels), and so the scheduling policy maintains fairness among users. 

Up to now, we have obtained the general form of the sum rate and individual user rate with the help 
of Gk( e ) without assuming specific distributions on the channel models. In the next section, we derive 
the closed form expression for Gk{e) with standard channel models. 
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IV. Exact Performance Analysis for Rayleigh Fading Channels 

In this section, we perform exact analysis to derive the closed form sum rate with standard Rayleigh 
fading channel models. Section HV-AI examines the PDF and CDF of the SINR for each user and uses 
them to derive a closed form expression for Qk{e) in Section HV-BI 

A. The Statistics of CQI 

In a practical system setting, the time scale for the large scale and small scale channel effects are much 
different. The variation of the small scale channel gain H occurs on the order of millisecond; whereas 
the large scale channel gain G which may consist of path loss, antenna gain, and shadowing, varies 
usually on the order of tens of seconds. Therefore, the large scale channel effect is assumed to be known 
in advance by the system, through infrequent feedback or location awareness. The small scale channel 
effect is modeled as complex Gaussian distributed random variables with zero mean and unit variance 
£A/"(0, 1). From the definition of SINR in ©, it can be seen that the numerator is a scaled x 2 {2) random 
variable (i.e., chi-square random variable with 2 degrees of freedom), and the denominator is a weighted 
sum of x 2 (2) random variables plus a constant. The following lemma provides the density function of 
SINR[ 0) , namely f z m. 



Lemma 1. The PDF of can be expressed as 



/ z (o,(x) =J2^ { k e 



1 i , gf/f x 



u(x), (15) 



m (fa) 

where w k = \\ {b p ) k <,> . and u(-) is the Heaviside step function. 



Proof: The proof is given in Appendix |A] 

7 jf\ namely F z{ . 



From Lemma [B the CDF of zf \ namely F„(o) can be computed as 



F zW (x)= / > w™e V m m + 

Z fe /n ^ k \ M l J 6 )„, / (0) . (b) 




i , (0) (6) 

1 , Pk Pk 



(b) (0) 
7 k e Pk Pk 

~1 Pk ] +Pk )x 



dy 

Pk' + Pk'v (pT+pVv)~) 

u{x). (16) 



B. Procedures to Compute Gk{t) 

Now we consider the computation of Gk( e ) = /o°° l°g2(l + x )^(^ ? z <0) wmc h Wli l be carried out in 
three steps. Step 1 provides a suitable PDF decomposition of d(F z (<» (x)) e by examining the expression 
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for the PDF. In Step 2, the decomposed PDF is further expanded for integration. Finally, Step 3 employs 
the outcome of Step 1 and 2 to derive the closed form expression for Qk{e) by standard integration 
techniques. The details are presented next. 

Step 1: We are interested in the exact formulation of d{F z ( 0) {x)) e , where the exponent e S N+. In the 
following lemma, an amenable decomposition is proposed for the statistical form of d(F z m (x)) e . 

Lemma 2. The PDF d(F z ( Q ) (x)) € with e 6 N + can be decomposed as 



e-l 



d(F zi0) (x)y = eJ2 



e-l\ (-1)' 



/ 



+ 1 



d 



1 - e 



V 



1 



J fe (0) 



M k (M (0) 
6=1 P k X + % 

Pk 



\ 



£+1 



J 



(17) 



Proof: The proof is given in Appendix |A] ■ 
Step 2: Even though the complicated form of d(F z ( )(x)) e is decomposed into (fTTT ). its formulation 
still prevents direct integration. The following lemma provides an expanded form for one of the terms to 
facilitate further integration. 



Lemma 3. 



J k (o) 



1 



6=1 




E 



3i+-+3J k 



Ji,---,3J k 



(b) ( 

\ j* * V'fc,- n ( 



= 1 4 = 



where 



0, i = 

1 d j b~ 

(jb— i)! dx 3 t>- 



0) \ jb Jfc 



6=1 



(o) \ -> 



Pi, 



(0) \ « 

Pic 



t > 1. 



(18) 



(19) 



"F 7 



Proof: The proof is given in Appendix |A] ■ 
For illustration purpose, the formulation of Lemma |3]is discussed and provided for the special Jk = 2 
case in Appendix lAl 

Step 3: The following theorem completes the final step by utilizing the outcomes of the above two 
steps to derive the closed form expression for £7&(e). 

Theorem 1. Gk( e ) can be computed as 



e-l 



£-l\(-l)< 



+ i 



E 



n 

6=1 



h+~+3J k =e+i 
C' + i) pI 1 



£ + 1 

31, ■ ■ ■ Jj k 



Jk jb 



EE< 

6=1 i=0 



(b) 
% 



(0) 



(0) ' (6) 
Pk Pk 



(20) 
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where T\(a, (3, 7) = J °° ^ x wnose closed form expression is presented in Appendix lAl 

Proof: The proof is given in Appendix |A] ■ 
The three-step procedure yields the closed form expression for £7fc(e), which can be substituted into 
(fT2l and ([141 to compute the closed form sum rate and individual user rate. The exact closed form 
expressions only involves finite sums and factorials making it computationally tractable and useful for 
system evaluation. In the following, the treatment of two simplified special cases are provided: the one- 
dominant interference limited case and the noise limited case. 

One-Dominant Interference Limited Case: This case approximates the scenario when there is one 
dominant interferer. Without loss of generality, assume p^' 3> for b ^ 1 and the effect of noise is 

(Q) (0) f (0) \H m 2 

omitted. Then the SINR can be approximated as SINR^ ~ SIR^ = pjTOTTii , which is the F distributed 



random variable. The CDF of the CQI can be written as 



F z p(x)=(l- qA- )u(x). (21) 
In this case, the computation of Gk(e) can be reduced to 

0- W j^f E ('-, ') ^Be,a(U + (l. l;l + 2; 1 - J) , ,22) 

where (a) follows from (371 3.197.5], Beb&(x,y) = /q 1 ** -1 ^ - J)" -1 * is the Beta function (57j 8.38] 
and 2-^1 (-, •; •; •) is the Gaussian hypergeometric function ||58l . 

Noise Limited Case: This case approximates the scenario when the impact of intercell interference is 
negligible, i.e., < 1. The SINR can be approximated as SINR*™ ~ SNR^ = pf^H^J 2 , which is 
the x 2 (2) distributed random variable. The CDF of the CQI can be written as 



F 7m (x) = 1-e p * u(x). (23) 




In this case, the computation of Qk{e) can be reduced to 

where (a) follows from 11571 4.337.2] and E\{x) = ^j-dt is the exponential integral function of the 
first order ||58l . The noise limited case is equivalent to the single cell problem which has been addressed 
extensively in HD, EH, 11461 . 

Remark: It can be easily seen that the CDF-based scheduling for the two simplified cases has the 

JO) 

same effect as the "normalized" CQI based scheduling, which is normalized by ^fa for the one-dominant 

Pk 

interference limited case and normalized by pl°' for the noise limited case. The general CDF-based 
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scheduling policy enables the general analysis for the multicell scenario, whose closed form expressions 
have been obtained by the aforementioned procedures. The exact expression, though computable, is not 
easy to interpret and draw insights. We now use asymptotic analysis to develop results that have the 
potential of providing further insights. 

V. Asymptotic Performance Analysis 

This section is devoted to the asymptotic analysis when the associated users in a given cell grows 
large. Section IV- A I proves the type of convergence exhibited by the received CQI under best-M partial 
feedback. A brief summary on the different types of convergence is provided in Appendix |B] for easy 
reference. In Section IV-B I the asymptotic rate approximation is derived and is employed to determine the 
minimum required partial feedback in Section IV-CI The results are presented with the proofs relegated 
to the appendix. 

A. The Type of Convergence 

The first step towards performing asymptotic analysis is examining the tail behavior of the received 

CQI at the scheduler side under partial feedback, namely for user k. In the full feedback case, 

F Y (o) = F z ( ) , which means the tail behavior of the CQI at the user side is equivalent to that of the 

5^ jy-. However, for the general best-M partial feedback, the relationship between Y^°j^ and is given 

by ® and is recalled here for easy reference as F y(0 , (x) = Y^Zl Si(N, M, m){F z(m {x)) N - m . One 

natural question is concerning the relationship between the tail behavior of Y^L and or formulated 

in a rigorous way: how to infer the type of convergence of F v m from the type of convergence of F 7 { ) 

r fc,M " " bl- 

under the condition of best-M partial feedback? The following theorem addresses this issue. 

Theorem 2. (Type of Convergence under Partial Feedback) F v (o) has the same type of convergence 
property as i^c) under the best-M partial feedback strategy. 

Proof: The proof is given in Appendix |B] ■ 
Theorem [2] states that the best-M partial feedback does not affect the type of convergence. In other 
words, once the type of convergence for F„< ) is proven, the same property is established for F v ( ) . Note 

' *k,M 

that so far no specific statistical property has been assumed for F z ( ). In the following, the statistical 
model expressed in (fl6l ) will be utilized for further analysis. The following corollary describes the tail 
behavior of F 7 <o> . 
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Corollary 1, For the general SINR case in the multicell scenario, F 7 m and F v w belong to the domain 

^k *k,M 

of attraction of the Gumbel distribution iPTTl . i.e., F 7 ( ) G P(G3) and F v (o) G "D(Gs). 

^/s *k,M 

Proof: The proof is given in Appendix |B] ■ 
For completeness, the tail behavior of F z (o) for the special simplified one-dominant interference limited 
case and the noise limited case is also provided in the following corollary. 

Corollary 2. For the one-dominant interference limited case, F 7 (o) and F v ( ) belong to the domain of 
attraction of the Frechet distribution P7l . i.e., F~m G T>{G\) and F Y (o) G "D(Gi). For the noise limited 

^k *k,M 

case, F z ( ) and i^y(o) belong to the domain of attraction of Gumbel distribution ll47l . i.e., F„( ) G T>(Gs) 

k k , M k 

andF v ( , G X>(6? 3 ). 

Proof: The proof is given in Appendix |B] ■ 
These established type of convergence results will be used to obtain the asymptotic rate approximation. 



B. Asymptotic Rate Approximation 

We now investigate the asymptotic approximation for the exact sum rate whose closed form expression 
has been derived in Section JV] Two additional issues arise in the heterogeneous multicell setting under 
partial feedback when compared with the standard homogeneous setting under full feedback. 

The first issue regards the heterogeneous statistics of the SINR for different users. In the homogeneous 
setting, the maximization or the order statistics is over the same CDF. Recall that the use of CDF- 
based scheduling in this paper has enabled each user to virtually feel that the other associated users are 
experiencing the same CDF for scheduling competition. Therefore, for a given user k, the order statistics 
is over the CDF of user fc's received CQI, which makes the individual user rate more interesting than 
the sum rate. 

The second issue arises due to the effect of partial feedback. In the full feedback case, the number 
of CQI values to maximize over at the scheduler is fixed and equals the number of the associated users 
K. However, due to partial feedback, the number of CQI values to maximize over at the scheduler is a 
random quantity. In other words, partial feedback results in a random effect on multiuser diversity. In the 
exact analysis in Section Jill this effect is reflected in the use of U n> M- We are interested in the asymptotic 
effect when the number of users grows large. To examine this random effect in the asymptotic analysis, 
denote the sequence of random variables K n {K) as the number of CQI values fed back for resource 
block n with K associated users. It is easy to see from © and © that n n (K) are binomial distributed 
with probability of success 4r under best-M partial feedback. Thus by the strong law of large numbers, 
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as K grows, the number of CQI values fed back for each resource block becomes Moreover, the 
convergence property of the sequence K n (K) can be shown by invoking the central limit theorem 



where d indicates convergence in distribution. Therefore, by employing the techniques which study the 
extremes over random sample size iBTI . ll60l . we have the following lemma. 

Lemma 4. When the number of associated users K goes large, the extreme order statistics lETTI of the 

KM 

received CQI for a given user k can be efficiently approximated by ( F v ( ) ) N . 

Proof: The proof is given in Appendix |B] ■ 
Now consider the limiting distribution of the maximum rate in order to derive the asymptotic approx- 
imation for the exact rate. Specially, we examine the limiting distribution of the rate i?fc,M> 

Rk,M = T(Y k}M ) = log 2 (l + Y k)M ), (26) 

where the superscript (0) is temporally dropped for representation simplicity, and will be added later to 
tailor the results for specific Kq and yS },. Note that the function T(-) in (1261 ) makes it tedious to directly 
check the conditions needed to enable finding the form of the asymptotic distribution. In |[50l . a limiting 
throughput distribution theorem is proposed for the full feedback single cell case for a narrowband system. 
Herein, we generalize the result to be applicable to the general SINR case in the general partial feedback 
OFDMA scenario with the following best-M limiting throughput distribution (LTD-M) theorem. 

Theorem 3. {LTD-M Theorem) Assume that under the best-M partial feedback strategy with N resource 
blocks and K associated users, the CQI received at the scheduler for user k, Y^^m is a nonnegative 
random variable with CDF Fy kM (x) such that fy k M (x) = F y (x) > and uj(Fy k M ) = sup{x : 
F Yk M (x) < 1} = oo. If lim j^^. = <f> > 0, F Yk M € V(Gi), i.e., F Yk M belongs to the domain of 



X— )-00 x - 1 Y k,M 1 



attraction of the Frechet distribution, or if lim 4- 



-jz — nr - 0, F Y e V(G 3 ), i.e., Fy belongs 
to the domain of attraction of the Gumbel distribution, then the distribution of the throughput for user 
k, Fji k M (r) = Fy k M (T _1 (r)) G V(Gs), i.e., Fn h M belongs to the domain of attraction of the Gumbel 
distribution. Moreover, the normalizing constants P71 for user k are given by 

a k:K (M) =lo : ' ' "-' ' ' A 

b k:K (M) = log 2 ( - ; ^ \ — ^ ) . (27) 
Proof: The proof is given in Appendix |B] 
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Remark: The LTD-M theorem enables us to study the distribution of Yfc,M instead of directly examining 
Fr u m . Also, note that the relationship of the type of convergence between Y^a/ and has been revealed 
in Theorem [2] Thus the connection between and Rk,M can be established by combining the two 
theorems. 

The normalizing constants in Theorem [3] can be used to obtain the asymptotic rate approximation. 
Denote c[ (M) as the asymptotic approximation for the individual rate of user k in cell Bq with total 
associated users Kq. Then according to the property that convergence in distribution for the maximum 
nonnegative random variables results in moment convergence |[50l , floTTl , C^{M) can be evaluated by 
the normalizing constants as follow^] 

cf\M) - ± (l - (l - f)"") („™ S M) + BMgjiO) • (28) 

where Eq is the Euler constant, and (1 — M-) K is the probability of scheduling outage. According to the 
CDF-based scheduling policy, the asymptotic approximation for the sum rate, denoted by C^°\M), can 
be computed as 

C (0) (M) = ± (l- (l- f )*°) p (41(M) +^e W) ■ (29) 

The form in (1291 is simpler than the exact analytic expression derived in Section JV] and can be an 
alternate basis for studying heterogeneous networks. Looking again at the normalizing constants in (l27l ). 
the specific expressions involve the inverse of the distribution function, Fy (•). In general, due to the 
complicated form of the SINR as well as the procedure to evaluate Fy k M , this inverse function can not be 
expressed in simple closed form except in some simplified cases. Since the CDF is a function of a scalar 
and is monotonically increasing, standard iterative algorithms are well suited for its computation. Now 
we consider the two aforementioned simplified cases: the one -dominant interference limited case and the 
noise limited case for illustration. For these two special cases with full feedback and best-1 feedback, 
the inverse of the distribution function can be computed in closed form, which are summarized in the 
following corollaries. 

Corollary 3. In the one-dominant interference limited case under full feedback, the specific form for the 

9 This form of asymptotic approximation leverages the first and second order moments of the extreme order statistics. Dealing 
with higher order moments and eventually the rate of convergence is left to our future work. 
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normalizing constants are given by: 



a 



(o) 



(0) 

Pk 



i:Ko (N)=log 2 ^1 + ^(^-1)1, 



41(^0= log 



Pk 

(0) 

1 + ^(^0-1) 



(30) 



In the noise limited case under full feedback, the specific form for the normalizing constants are given 
by: 



(0) 

Pk 



l + pf lnif 0/ 



(31) 



Proof: The proof is given in Appendix |B] ■ 

Corollary 4. In the one-dominant interference limited case under best-1 feedback, the specific form for 
the normalizing constants are given by: 

(K -N)^ \ 



41(1) =^211 + % „ 

V Pi' K» -(Kq-N)s 



(0) 

Pk 



4°k(l) = l°g 2 



(l 



\ 



' (1 



(K e-N)Tr 
— t 



Pk {K Q e)N-(K Q e-N)N 



1 -I- ?k 



(Ko-N)Tr 



(32) 



In the noise limited case under best- 1 feedback, the specific form for the normalizing constants are given 
by: 



,(o) 



^ ( 1 )= 1 °g 2 ( 1 + 4° ) ln — 



6 (0) m 



log 2 



K» -(K -N)- 

(K e)TT 
(K e)7T-(K e-N)w 



\ 



K' 



(33) 



J 



Proof: The proof is given in Appendix |B] ■ 
For the general best-M partial feedback case, the normalizing constants can be obtained using (T27T ) and 
the specific CDF for the corresponding case. To illustrate the benefit of asymptotic analysis for the two 
simplified cases, we conduct a numerical study to compare the sum rate obtained using the exact analysis 
and the asymptotic one under different symmetric large scale effects in Fig. [2] It is interesting to note 
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Number of Associated Users Number of Associated Users 

(a) (b) 

Fig. 2. Comparison of the sum rate for best-M feedback obtained using the exact analysis and the asymptotic analysis under 
different symmetric large scale effects for different M with respect to the number of users (TV = 16; M — 1,2,4): (a) the 
one-dominant interference limited case (p = ^nr); 0°) me noise limited case (p = p' '). 



TABLE I 

The Main Steps of the Analytical Framework for the Exact Analysis and the Asymptotic Analysis 



Framework 


Exact Analysis 


Asymptotic Analysis 


Technical Steps 


Main Results 


Technical Steps 


Main Results 


Step 1 


Statistics of zf ] 


Lemma 1 — 3 
Theorem 1 


Type of Convergence of F (o) 


Corollary 1 


Step 2 


Statistics of Y^ M and y£ 0) m 


Type of Convergence of F ' ( ) 

'k.M 


Theorem 2 


Step 3 


Statistics of X$ \ k* = k, \U M \ = r 


Type of Convergence of F (o) 


Theorem 3 


Step 4 


C (0) (M) and cf } (M) 


C ((,) (M) and C ( °\M) 


Theorem 3 



that the asymptotic expressions hold well even for small number of users, which means the convergence 
to the limiting distribution is fast. 

Up to now, we have leveraged exact analysis and asymptotic analysis to derive useful closed form 
results for the exact sum rate and the asymptotic approximation. For easy reference, the main technical 
steps and the main results are summarized in Table U In the next part, the procedure to determine the 
minimum required partial feedback is examined. 

C. Determining the Minimum Required Partial Feedback 

Firstly, the asymptotic optimality of the best-1 feedback is presented in the following corollary. 
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Corollary 5. When the number of associated users K — > oo, the performance loss of using best-1 feedback 
in terms of sum rate approaches zero. 

Proof: The proof is given in Appendix |B] ■ 
We are more interested in the pre-asymptotic user regime where the simple best-1 strategy is no longer 
optimal. The goal is to choose the minimum required M without seriously degrading system sum rate 
when compared to a system with full feedback. The selection of M can be formulated as the solution to 
the following optimization problem: 

C(°)(M*) 

Find the minimum M , s.t. — 7— — — > 77, (34) 

C(°)(JV) 

where (M) refers to the exact analytic expression derived in Section [TV] and rj is a system defined 
threshold. As mentioned previously, the established asymptotic expressions are more computationally 
efficient than the exact one. By leveraging the asymptotic approximation, the problem (l34l can be 
reformulated as 

C(°)(M*) 

Find the minimum M* , s.t. — 7— — — > 77, (35) 

C(°)(iV) 

where C^(M) refers to the asymptotic approximation in (|29l l. 

We see from (l34l and (I35T ) that M* or M* depends on the number of users associated with the base 
station and the corresponding large scale channel effects. Since these factors can be highly diverse in 
a heterogeneous network, the minimum required partial feedback is inherently different across different 
cells. Therefore, by employing the established asymptotic results, M* can be quickly determined to track 
M* in order to design situational-aware heterogeneous partial feedback. 

VI. Numerical Results 

In this section, we conduct a numerical study to support our analysis. A simple heterogeneous network 
is modeled with two macrocells each with two picocells inside. The locations of the picocells are randomly 
placed and then fixed for simulation. The system bandwidth is 5 MHz, the noise power spectral density 
is —170 dBm/Hz, and the number of resource blocks N = 16. The transmit powers of the macrocell 
and picocell are 43 dBm and 30 dBm respectively. The path loss (in dB) model in [62] with 2 GHz 
central frequency is employed: the path loss from the macrocell base station to users is 15.3 + 37.6 log 10 d 
for distance d in meters; the path loss from the picocell base station to users is 30.6 + 36.71og 10 d for 
distance d in meters. Log-normal shadowing is assumed with standard deviation of 8 dB. The radius of 
the macrocell and picocell is assumed to be 500 m and 100 m respectively. For each drop in simulation, 
users are randomly placed and the cell association is determined by the large scale effects and fixed. 
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— H — Greedy Scheduling 

O 1 1 Round Robin Scheduling 
-A- CDF-based Scheduling 
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— B — Greedy Scheduling 
O Round Robin Scheduling 
A CDF-based Scheduling 
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(a) 



(b) 



Fig. 3. Comparison of the three scheduling policies: the greedy policy, the round robin policy, and the CDF-based policy 
under best-M partial feedback strategy with respect to the associated users in a macrocell (N — 16; M — 4): (a) The sum rate 
comparison; (b) The system fairness comparison. 



Firstly, the performance of CDF-based scheduling is compared with the greedy scheduling policy and 
the round robin scheduling policy in Fig. [3] In this simulation, users are assumed to employ the best- 
M partial feedback with M = 4. In Fig. [3] (a), the sum rate with respect to the number of associated 
users in a macrocell is shown, which is averaged by performing 1000 independent drops. It can be seen 
that the round robin policy does not invoke multiuser diversity at all, and the sum rate performance 
of the CDF-based policy is close to the greedy policy. Fig. [3] (b) compares the system fairness for the 
three policies. The system fairness is defined and discussed in |[63l , |[64l using the following form: 

= — 'ff , where P& refers to the proportion of resources assigned to user k with the 
normalization factor In if. The system fairness for the round robin policy and the CDF-based policy is 

1 despite the number of associated users and the heterogeneous channel effects. However, the system 
fairness for the greedy policy is decreasing when more users are associated. This is due to the fact that 
some high geometry users occupy the system resources with high probability when the greedy system 
has more serving users. Therefore, the CDF-based scheduling policy enjoys the multiuser diversity while 
guaranteeing fairness at the same time, which makes it well suited for the heterogeneous framework. 

Next, in order to evaluate the derived closed form results from exact analysis and the corresponding 
asymptotic approximation, one individual user is randomly selected for demonstration. To illustrate the 
scaling performance with respect to the number of associated users, the so called individual sum rate 



TO APPEAR IN IEEE TRANSACTIONS ON SIGNAL PROCESSING 



21 



7 




Fig. 4. Comparison of the individual sum rate (individual user rate multiplied by the number of associated users) obtained 
from simulation, analytical result, and asymptotic approximation for different partial feedback M with respect to the number of 
associated users for a randomly selected user in macrocell (TV = 16; M = 2, 4, 16). 



for this user is of interest. The individual sum rate is the individual user rate multiplied by the number 
of associated users. Fig. @] plots the individual sum rate obtained from the analytical expression by 
exact analysis, from simulation, and from the established asymptotic approximation using normalizing 
constants. Different partial feedback cases are also shown for comparison. It can be observed that the 
analytical expression and the simulation results are in full agreement. Also, the asymptotic approximation 
tracks the system performance very well. Furthermore, the rate gap between the partial feedback case 
and full feedback case becomes negligible when the number of associated users increases. 

Finally, Fig. [5] (a) examines the minimum required partial feedback M* obtained by using the expression 
for the sum rate from the exact analysis and the asymptotic approximation for a macrocell. Two different 
thresholds are set for evaluation: r] = 0.9, 0.99. It can be seen that the results obtained using asymptotic 
analysis track the results from exact analysis very well, especially for lower threshold and larger number 
of users. Since the number of users associated with each cell as well as the large scale channel effects can 
lie in diverse ranges, this results in heterogeneous partial feedback in heterogeneous multicell networks. 
The total number of partial feedback with respect to the number of associated users for a macrocell is 
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- Asymptotic Analysis (n=0.9) 

- Exact Analysis (t)=0.9) 
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Fig. 5. Partial feedback with respect to the number of associated users for a macrocell (N — 16): (a) Comparison of the minimum 
required M obtained from exact analysis and asymptotic analysis under different thresholds (77 = 0.9, 0.99); (b) Comparison of 
the total number of required partial feedback from exact analysis and asymptotic analysis under threshold r\ = 0.9. 



illustrated in Fig. [5] (b) under the threshold rj = 0.9. The total number of partial feedback is calculated 
by multiplying the minimum required partial feedback M* with a given number of associated users. It 
can be seen that the range of variation for the total number of partial feedback is limited. Even though 
the number of associated users is 5 times larger, the total number of partial feedback does not change 
much. This is due to the heterogeneous partial feedback design to adapt the number of partial feedback 
to the number of users as well as the channel conditions. 

VII. Conclusion 

In this paper, an analytical framework is proposed and developed to investigate the performance of 
situational-aware heterogeneous partial feedback in an OFDMA-based heterogeneous multicell using the 
best-M partial feedback strategy. The system model is general and thus the obtained results can be 
generalized and applied to conduct system evaluation with alternate statistical models. The CDF-based 
scheduling policy employed in this paper has the desired property of supporting multiuser diversity 
while maintaining scheduling fairness among the contending users to guarantee each user's data rate 
despite of different locations and large scale channel effects. The exact closed form sum rate is obtained 
for the multicell model by suitable decomposition and expansion of the received CQI at the scheduler 
side. Asymptotic analysis is carried out to draw further insight into the multicell model with partial 
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feedback. Interestingly, the effect of partial feedback does not alter the type of convergence of the 
received CQI. The random effect of multiuser diversity caused by partial feedback is also examined and 
asymptotic approximations are derived by utilizing the normalizing constants. The established asymptotic 
approximation tracks the exact system performance well even for small number of users. Therefore, it 
can be leveraged to quickly determine the minimum required partial feedback in a given cell. In this 
paper, only a single antenna is assumed and the scheduling policy is fully distributed. Our future work 
will consider a multi-antenna setting and more advanced scheduling policies. 



Appendix A 

Proof of Lemma\l} Denote £ = Y2b=i Ph\ H kl\ 2 ^ which is a we ig hte d sum of x 2 (2) random variable. 
In practical system setting, the large scale effects from the interfering cells are distinct, the PDF of ( is 
derived to be [65], 11661 : 



where w 



(b) 



W Pk J 

(b) 

II (bf k (i) ■ Then the PDF of Z k can be obtained as follows 

i=i Pk ^Pk 

i^b 



f m (x) = f z ( 0) ,Jx\()f c (C)d( 

k J Q fc is 

i » J h „( 6 ) poo 

Pk 6=1 Pk J0 



x Jk (6) 

„(0) e (6) 
Pfc 6=1 /'fc 



/ 



(0) 
Pk 



+ 



(0) 

Pk 



(0) \ 2 

A' 
(i. 
Pk 



u(x). 



(36) 



(37) 



Proof of Lemma EJ- It is clear that d(F z (o) (x)) e = e(F z (a)(x)Y 1 f z (o)(x)dx. Employing the binomial 
theorem [67], and substituting the expressions for F z m(x) and f z ^{x) yield 



e-1 / lX / .A (6) ~7W (0) 

' \ 6=1 Pfc + Pk X 



( 



1=0 



Jk (b) n 



Jo) Jb)(b)-^m 
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+ 1 



/ 



V 



Jk „{b)-^) n {0)\ \ 



— (0) . (6) 
6=1 Pk +Pk X 



(38) 



where (a) follows from the differentiation property of (1 — F z <.o)(x)) +l . 
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Proof of Lemma \3\ Applying the multinomial theorem ll67l yields 



J k (o) 



e+i 



1 



6=1 



(*>) o m 

Pk 



E 



e + i 
n, ■ ■ ■ Jj k 



n 



P i b} 



(39) 



Exploiting the partial fraction expansion [57 1 generates the expanded form in dl~8T > with specific expression 
for ifjf] defined in ([T9l . 

For illustration purpose, the expansion for Jk = 2 case, which corresponds to two major interferers is 
now presented. In this case, applying binomial theorem is sufficient for expansion, which yields 
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trivial i > case. 
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Proof of Theorem\J} The outcomes of Step 1 and Step 2 lead to direct integration to calculate Gk(t) 
as follows 
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where (a) follows from integration by parts. The form of (l20l is expressed by the definition Ii(a, (3, 7) = 
Jo°° (i+x)(p+x)~t d x - in order to compute Ti(a, j3, 7) into closed form, firstly apply partial fraction expan- 
sion 1571 to ^ l+x ^ +x y, , and then define the auxiliary integration 12(0, P, 7) — / °° (p+x)-< ^ * s c l ear 
that when 7 = 0, Zi(cn, /3,7) = ^(a, 1, 1). For the non-trivial case when 7 > 1, employing the partial 
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fraction expansion and after some manipulation the following relationship is revealed between I±(-, •, •) 
and X 2 (-, -) : 

Ii («, ft 7) = (j g^ 1)7 J 2(a, 1, 1) + E {r^ X2(a ' ^ 7 " 1 + 1} ' (42) 

Z 2 (a,/3,7) can be further computed by noting from ll57l 3.352.2] that X 2 (a, /?, 1) = e a/3 E\(af3), where 
E\{x) = J°° ^j-dt is the exponential integral function of the first order ll58l . and utilizing integration 
by parts. The closed form result for X 2 (a,/3,7) is presented as follows 

( _l)T-i^-i e ^ El ( Q/ 3) (,_!)! s _! _! 

X 2 (a,/3, 7 ) = < (7 lj - *=i (7 1} - (43) 

e a/3 Fi(a/3), 7 = 1 

Appendix B 

Lemma 5. (Sufficient Conditions for Type of Convergence / I47I/ - 0O1/ ) Let Ai, A 2 , • • • , A^ be i.i.d. random 
variables with CDF F\(x). We denote A# = maxj Aj. If there exists some distribution function G which 
is nondegenerate and some constant cik E M, 6 a' > such that the distribution of AK bK aK converges to G, 
then G must be one of just three types: G\. Frechet distribution; G 2 : Weibull distribution; Gy. Gumbel 
distribution. 

The CDF of Aj, i.e., F\ determines one of the exact types. If F\ results in one limiting distribution, then 
we say Fx belongs to the domain of attraction of this type, i.e., F\ G V(Gi). The well-known sufficient 
conditions for Fx G T>(G\) and Fx G T>(Gs) are as follows: Define to (Fx) = sup{x : Fx(x) < 1}. F\(x) 
is absolutely continuous and f\(x) = F^(x) and f'\(x) = F"(x) exist, then 

(a) Fx G V(G\) if f\(x) > for all large x and for some <fi > 0, 

x/a(x) = 0> 

z->oo 1 — fx(X) 

(b) Fx G T>(G-2) if /u < oo and for some > 0, 

(fj,-x)f x (x) 

lim — ? r = 0. (45) 

1 - Fx(x) 

(c) Fx G T>(G?j) if /a(x) > and is differentiable for all x in (xi,oj(Fx)) for some x\, and 



lim — — 

x^u(Fx) ax 



0. (46) 



fx(x) 

Further, we can choose the normalizing constants ax = F\ (1~ -4) and 6^- = F x ~ 1 (l — ^-)— FT (1— -4), 
where F x _1 (x) = inf{y : F\(y) > x}. 
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Proof of Theorem |2} Assume zj® is a nonnegative random variable with CDF F z ( ) (x) such that 
f z (o) (x) > and /^ (0) (x) exist. The random variable Y^ ^ is related to zj® by the following equation: 



^M-l 



F v ( 0) (x) = YjZ- r \ii{N,M,m)(F 7 <. )(x)) ■ in order to snow that F vi0) has the same type of 

*k,M *k,M 

convergence property as F z (a) , the proof in the sequel will be conducted for each of the three types. 

s=/ z ( )(aO x f Y (S>) 0*0 

(i) If for some 6 > 0, lim , — =* — ^ = 6, then F 7 (oj G It must be shown that lim , — „ k ' M , , = 

for some 6 > 0. Substituting the expression for F v <o> and / v ( ) yields 



lim 

x— >oo 1 



x/ y( o) (x) 



* Em=o eiC^V, M m)(iV - m)(F z ( 0) (x))^— (x) 



lim 



> , « 0*0 ™ 1 - E™=o 6(iV, M, m)(F z <o, (x))^-™ 



(a) 



lim 



Em=o Zi(N, M, m)(N - m) (/>) (x) + xf m (x) 



^-E^=oei(^M,m)m' 
where (a) holds by applying the L'Hospital's rule; (b) follows from the type of convergence of Z£' '. 
Therefore, 6 = 6, and F v < ) G T>(G\). 

(ii) If u < oo and for some > 0, lim -,_ F — fc , t = 0, then F 7 (o) G DOG?). It must be shown that 

X— >/J 1 P Z («)\ X ) ^k 

(P-X)fy(.0) (X) 

nm i_f — fc ' A L\ = f or some > 0. Substituting the expression for F v (o) and / v (o> yields 



(47) 

(0) 



(/i - x)/ y( o) (x) 



- j E^=o MjV, M, m)(iV - m)(F z (o, (x))^- 1 /^, (x) 
1 " E™=o £i W M, m)(F z ( 0) (x))^— 
E„=o 6(iV, M, m)(N - m)(F z m (x)) w — 1 (/^ 0) (x) - ( M - x)/;< 0) (x)) 



lim 

1 — is 



1 i. 



(x) 



lim 

X— >/i 



= lim ■ 

x— >n 



E™=o M, m)(N - m)(F z w (x)) 



N—m—l 



/ Z (0) (x) 



(6) E™=o CiW M, m)(7V - m)(i> 5 (x)) 
" lim 



N-m-l 



(48) 



c ^ Em=o M > ™)(N ~ m)(F z m (x))^-™- 1 
where (a) holds by considering the term that dominant the limit and applying the L'Hospital's rule; (b) 



follows from the type of convergence of zP . Therefore, 



, and F Y w G V(G 2 ). 



(hi) If lim A 



lim 



tlx 



1-F 



•(0) 
k.M 



1-F (0) (x) 
71(0) (a:) 



1-F z<0 )(x) 

0, and lim — > — V^— exists, then F„(oj G V(G^). It must be shown that 



/ v (o) (s) 

fc,M 



0. Carrying out the differentiation, another equivalent condition is the following: 
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K (o) 

k . M 



/ k.M 



i^o (/ (o) (a:)) 



1. Substituting the expression for F Y (o) and f Y ( ) yields 



lim 

x— >oo 



v (o> (x)-l\f (x) 

k,M / Y k.M 



lim 



( E^o 1 6(Ar,M,m)(F z(0) (x))»-™-l) fEl^o 1 6 (iV,M, m )(JV-m)(iV-m-l)(F (0 , (x))"-»- 2 (/ z < 0) (x)) : 

\ fc / \ k k 



+ lim 

x— >oo 



I Em=o e 1 (iV,M,m)(iV-m)(F z (o ) (x))*-"- 1 /^) 

\ k k J 

(Ef^o 1 ttW^f) (E^o 1 ei(JV,M,m)(/V-m)(F z ( 0) (x))« — V; m (*)) 



W (EL-o'^(W(^)); (E:::,'a^^)(f r ( I ))^"-i)(ES6(^M, m )(N- m )(F r , M ) 

~~ ^m-i 7v fl/f m \nv—™\\' 2 (/ (0 )(x)) 2 



(6) 



(E^o6(^,M,m)(AT- m )) 



where (a) comes from the fact that ^m^ 1 £i(-^> m ) = 1> (b) holds by applying the L'Hospital's rule 
and the type of convergence of Zy. Therefore, F v <a) E T>{G^). 



(49) 



Proof of Corollary [7} For the general SINR case, in order to prove that F„(o) E P(Gs), it must be 

0. Substituting the expression for f z m and F^coj in (1151 ) and (1161 ) 



shown that lim 



yields 



/_(0) (*) 



lim ■ 

x— >oo 



1 - F (o) (x) 
/ Z (0) (x) 



, (6) P ' (0) 



lim 



6=1 



pi 0) +pi b) - 



(6) P 

2^6=i ro fc e 



+ 



/5 <0, o (6) 



M T J k (b) i 
(a) ^ ^6=1 ^fc p W 

2^6=1 w fc 



(0) 
Pk > 



(50) 



where (a) follows from applying the L'Hospital's rule. It can be shown that 



/U(*)=E- 



6=1 



(0) . (6) 



(0) , (6) 
Pfc +Pk X 



+ 



(0) , (&) 
Pk +Pk X 



Utilizing the same technique as in (1501 . it can be easily shown that 

f z (0) ( x ) 1 



lim 

x— >oo 



(51) 



(52) 



(V CO) (*)-l)/' (0) (*) 

V k / k 



By combining the results of (1501 ) and (1521 ). the following equation holds lim r . — ^ 

x->oo U z (o) 

— rrr - 1, which proves the type of convergence of F„(o>. Applying Theorem |2] yields F v ( ) E V(G^). 

~Pk " k " Y k,M 
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Proof of Corollary |2]- In the one-dominant interference limited case, it is easy to verify that 



lim 



xf z m (x) 

k 

F z m (x) 



lim 



(1) 

Pk x 



> , (0) 
Pk X + Pk 



1. 



(53) 



x— >oo 1 

Thus, F 7 ( ) G D(Gi). Applying Theorem |2] yields F v (o) G 

In the noise limited case, it is easy to verify that F z m G X^G^), e.g., see [50]. Applying Theorem [2] 
yields F y(0) G V(G 3 ). 

1 k,M 

Proof of Lemma |?} To prove this lemma, the following theorem which discusses the extremes over 
random sample size is called upon. 



Theorem 4. {Random Observations Theorem / l?7l/ . / I60I/ ) Let, as K — > oo, 

i? is a positive random variable. Assume that there are sequences ax G M, bx > such that — ^ 
converges weakly to a nondegenerate distribution function G. Then, as K — > oo, 



$ in probability, where 



/oo 
G y (x)dP(-& < y) 
-oo 



(54) 



— oo 
k(K) 



M 



. Thus from the above random 



From the analysis in Section IV-BI when K — > oo, — *g- —r jr. 
observations theorem, the extreme order statistics of the received CQI for a given user k can be efficiently 



approximated by F v ( 

Proof of Theorem \3} According to the condition of domain of attraction, it must be shown that 



lim — 

r->oo dr 



if lim xfYk - M 



% = cj> > or lim A 

W a;— >oo 



0. 



(55) 



0. 



It is derived in [501 that 



lim — 

r->oo dr 



lim 

r— >oo 



1 (l-^^W))/^^- 1 ^))' 



lim 

r— >oo 



(i-F y>iM (r- 1 (r)))(r- 1 ) ,, (r)) 
/n,M(r- 1 W)((T- 1 )'(r)) 2 



lim — — 



If lim 



I/.mM _ 1 



hm — 



(;r) 



(56) 



X— >00 ""J'k.M' 

(x) — >■ oo as x — y oo, if lim 



, then Fy kM G P(Gi). Using L'Hospital's rule, for a function 6{x) such as 



a;— >oo 



g(») 
x 



j, then lim = 4. This leads to lim 



r(r) 



0. 



If lim 4- 



x— >oo 



(/.)' 



fc.M 



0, then Fy k M G T>{G^). Similarly, applying L'Hospital's rule yields 



lim -f 

r— »oo ar 



0. 
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Up to now, the sufficient conditions have been proved. From the analysis in Section IV-Bl on the random 
effect of multiuser diversity due to partial feedback, the number of CQI values fed back for each resource 
block becomes with high probability. Additionally, note that 

F rL (*) = T (Fy^ M (*)) = log 2 (l + F^ m (*)) , (57) 
thus the normalizing constants (l27l i can be obtained. 



Proof of Corollary\3} In the one-dominant interference limited case with full feedback, the normalizing 
constants can be obtained by using F~^ 0) (x) 



(o) 

Pk x 



In the noise limited case with full feedback, the normalizing constants can be obtained by using 

^ 1 o)W=/>i 0) ln T ^. 

1 k,N 

Proof of Corollary^ In the one-dominant interference limited case with full feedback, the normalizing 
constants can be obtained by using F~} 0) (x) = ,/? fc x j_ and the number of CQI values fed back per 

Y k,N pI'{1-xn) 

resource block equaling 

In the noise limited case with full feedback, the normalizing constants can be obtained by using 
F~} m (x) = pf^ In — and the number of CQI values fed back per resource block equaling 

Y k,N 1-XN 

Proof of Corollary\5} When K — > oo, the probability of scheduling outage (1 — jt) k — > 0. Therefore, 
from (|27l it should be shown that lim — > 1. By applying the L'Hospital's rule, the following 

if— >oo (1 — jf) N 

equivalent equation holds: 

lim £ ^ r ; = 1. (58) 

A'^cx>l_(l_|)^ 
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